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Abstract 



We propose a new mechanism for late cosmological baryon asymmetry in 
models with first order electroweak phase transition. Lepton asymmetry arises 
through the decay of particles produced out of equilibrium in bubble collisions 
and is converted into baryon asymmetry by sphalerons. Supersymmetric mod- 
els with explicitly broken i?-parity may provide a suitable framework for the 
implementation of this mechanism. 



The realization that the baryon (B) and the lepton (L) violating quantum 
effects in the standard electroweak theory are efficient at high temperature [1] 
has sparkled a widespread interest on the issue of late cosmological B asymmetry 
production. The problem is that any matter-antimatter asymmetry created at 
some superheavy scale [2] can be easily wiped out by the 5-violating quantum 
effects, which are abundantly in equilibrium throughout all the history of the 
early Universe until its temperature drops down to the Fermi scale range [1]. If 
this is indeed the case, then we are faced with the vital problem of originating a 
new cosmological B asymmetry at the scale of the Electro Weak Phase Transition 
(EWPT). 

A key ingredient for the development of a sizeable late cosmological B asym- 
metry is that the EPTW be of first order. This seems necessary for the imple- 
mentation of the so-called out of equilibrium condition which, together with B 
and CP violation, constitutes one of the essential requirements to generate a net 
AB [3] . There already exist several studies of AB production at the electroweak 
scale in the Standard Model (SM) or in minimal extensions of it [4]. In SM the 
Z?-violating quantum effects play a twofold role [5]: on one hand they are respon- 
sible for the washing of any preexisting B asymmetry, but, on the other hand, 
they are invoked for the production of the new AB at the Fermi scale. In gen- 
eral, in order for these scenarios at the electroweak scale to work, the 5-violating 
quantum effects must be sufficiently suppressed after the accomplishment of the 
EWPT, or else they can wipe out the B asymmetry once again. This condition 
translates into an upper bound on the mass of the light neutral physical Higgs 
[5,6]. This limit is in contrast with or dangerously close to the present LEP lower 
bound on Higgs masses [7], depending on the details of the different analyses and 
whether the SM [5] or the Minimal Supersymmetric Standard Model (MSSM) 



are considered [6]. 

In this paper we focus on a different situation, where at the Fermi scale some 
new B- and/or L-violating phenomenon is effective in addition to the B- and 
L-violating quantum effects of the SM. In particular, we consider the case in in 
which new L-violating interactions are responsible for producing a cosmological 
AL at the EWPT, while the B- and L-violating quantum effects are subsequently 
called into play in order to transform this AL into the cosmological AB []. Notice 
that in our approach the requirement concerning the quantum effects just after 
the of the EWPT is opposite to that previously imposed in SM and MSSM: 
namely, they must be still efficient in order to convert AL into AB. Hence, the 
dangerous upper bound on the physical Higgs mass is removed in our case. We 
provide an illustration of these ideas in the context of supersymmetric (SUSY) 
models [9] where .R-parity is explicitly broken in the leptonic sector [10]. Unless 
the violation is extremely tiny, the L-violating transitions originated by i?-parity 
breaking are in equilibrium down to temperatures at which the EWPT occurs. 
We show that the .R-violating interactions can originate an adequate AL in 
theories with a first order EWPT and such AL can be efficiently converted into 
AB by L?-and L-violating quantum effects still operative at the end of the phase 
transition. 

Whether or not the EWPT in SM is of first order is still an open issue much 
debated in the present literature [11]. The new twist in this discussion is repre- 
sented by the inclusion of infrared effects which contribute non-negligibly to the 
definition of the order of EWPT. It seems that, in any case, the transition may 
remain first order even after their inclusion, although possibly weaker than what 

appeared in previous analyses [12]. Even a weak first order EWPT can constitute 

^aryogenesis via Leptogenesis was first considered in [8], however, in our model, the lepton 
asymmetry and the out of equilibrium condition are obtained in completely different manners. 



a suitable framework for our proposal. 

We consider the one-loop effective potential coming from the SUSY model 
whose Higgs sector contains, in addition to the usual two Higgs doublet superfields 
H® and H%, one or more singlet superfields iV's [13], with a typical coupling in 
the superpotential ^HiH^Ni. The theory simplifies in the limit in which the 
SUSY breaking mass scale rh is much larger than the weak scale. Indeed, only 
one combination of the two neutral Higgs bosons H® and H® remains light. For 
instance, in the case of one singlet N 

h = cos(3 ReH° + sin (3 ReH$ + O (-) ReiV, (1) 



x 



where 



and 



Re^ ) = v u Re(f/ 2 °) = v 2 , tan/3 = — , (2) 



Re(iV) = x, v = sjv\ + vl (3) 

Here we have assumed x ^$> v. 

Hence, at energy scales lower than rh, the theory approximately reduces to the 
standard model with only one Higgs doublet. The one-loop potential at finite 
temperature V (h,T) can be written as the sum of the classical potential, the 
Coleman- Weinberg quantum corrections [14] and a temperature-dependent part 
[15]. Calling To and v (To) the temperature at which the potential is flat at the 
origin and the vacuum expectation value of h at T = To, respectively, we obtain 

V -P- = f , (4) 

where a and At are the coefficients of the cubic and quartic term in h, respectively. 
At can be written as At = A — \K , where A is the coefficient of the h 4 term in 
SM and K denotes the correction coming from the supersymmetric contribution 



to the effective potential. Given our assumption that fh 3> mw, we can expand 
the mass squared of each SUSY particle in power of fh [6] 

m 2 = fh 2 + g t h 2 + O ( — ) . (5) 



777/ 



Then K is found to be [6] 



'-?^ff- « 6 » 



The sum in (6) extends over all bosons (+) and fermions (— ) present in the theory 

with number of degrees of freedom iVj. 

A non zero value of v (To) /To signals a first order phase transition. Here comes 

the abovementioned crucial difference between baryogenesis schemes where AB 

is generated before the end of the EWPT and our scenario. In the former case, 

in order to avoid the cosmic AB be washed out soon after the phase transition, a 

lower bound on v (T ) /T must be imposed, v (T ) /T ^1.3 [5,16]. From eq. (4) it 

is easy to realize that this bound implies an upper bound on the Higgs boson mass 

in SM where, roughly, m^ is proportional to A 1 / 2 . In our scenario the anomalous 

electroweak processes must be still effective soon after the accomplishment of 

the EWPT and hence, no upper bound on m^ has to be imposed. However, the 

presence of a singlet N is still relevant also in our treatment. It was recently 

observed [17] that in SM the EWPT may proceed via percolation of subcritical 

bubbles, originated by thermal fluctuation, unless the Higgs boson mass, and 

consequently A, are small enough, m/i< 100 GeV, according to the recent estimate 

of Dine et al. [11]. Thus, if one decreases the value of At, the first order phase 

transition is made easier. Taking the result for K from ref. [6] for the case of the 

SUSY model with two singlets N\ and N 2 , one can see that it is indeed possible 

to lower Xt considerably with respect to its value in SM 0, without decreasing 
2 In the presence of only one singlet N, the fermionic contribution increases At- We have 



the value of nih- In fact, one must be careful not to decrease At too much, or 
else v (To) /To becomes too large and the anomalous electroweak phenomena are 
no longer operative after the EWPT. Taking nih ~ iTi t m 100 GeV, a value of 
\t ~ 2 x 10 -2 generates a suitable phase transition for the baryogenesis scenario. 
We come now to a short description of the dynamics of bubble nucleation 
and collision in a first order phase transition. At high temperatures, T ^> Tq, the 
Universe is in the symmetric phase and there exists a unique minimum of V(h, T) 
at (h) = 0. As temperatures lowers, a second minimum appears and becomes 
degenerate with the minimum at (h) = at [18] 

where u is the coefficient of the h 2 T 2 term in the potential. 

The Universe decays from its false vacuum state at (h) = by bubble nu- 
cleation. Nucleation takes place at T NUC ~ T ~ 150 GeV, if Ay ~ 2 x 10~ 2 . 
The difference between the false and the true vacuum energy densities at T ~ To 

reads 

4« 2 w . . . 

PVAC « q^T T o\ (8) 

and is transformed into potential energy in the bubble walls whose energy density 
is (see Enqvist et al. [18]) 



2 3/2 



n 



■', 



At T < Tjsruc bubbles keep growing. If they expand by a factor 7 since T^uc, 
their kinetic energy becomes, roughly speaking, 0(7) larger than the potential 
energy [19]. The total energy of the bubble reads 

E TOT w 47ri? 2 ?77, (10) 



checked that the inclusion of the bosonic contribution may decrease At enough only if one 
performs a fine-tuning. 



where R is the radius of the expanding bubble. 

When bubbles collide the energy (10) can be released with two different 
mechanisms: i) direct particle production due to quantum effects [20] and ii) 
generation of coherent scalar waves [19,20] which decay into light particles (such 
a production appears to be less useful for our scenario). In both cases the result- 
ing particles are obviously produced with a distribution far from the equilibrium 
distribution and this represents the crucial ingredient in our scenario to imple- 
ment the out of equilibrium condition necessary to originate a net AL. 

Our aim now is to estimate the number of particles of each species generated 
through the abovementioned direct particle production mechanism. The mean 
energy of a particle produced in the bubble collision is of order of the inverse 
thickness of the wall, (E) ~ A"; 1 = A _1 7, where A m 6\/2\t/ (olTq) is the size 



of the wall when it is formed and we are assuming that 7 ph 1/v/l — v%,, v w being 
the wall velocity [19]. 

The number density of particles of species x produced in the collision is of 
order 

n X~ fxPVAC^j (11) 

where fa paramatrizes the fraction of particles of species \ produced in the direct 
bubble collision. We proceed to give an estimate of this factor fa. Following ref. 
[20], one computes the probability of producing particles x's in the decay of off- 
shell Higgs fields describing the dynamics of the bubble collision. Such probability 
reads 

v "i~Bf inm2 Hk2) e ^ 2 ■ 4m| ) ' (i2) 

where k = (k ,k x ,k y ,k z ), Ti(k) denotes the Fourier transform of the suitable 
Higgs field configuration, T{k 2 } is the imaginary part of the two points one- 
particle irreducible Green function constructed from the h-\-X coupling and m^ 

6 



represents the mass of the particle x- The important message of (12) is that 
particles x's with masses larger than rrih/2 can be produced through the decays 
of virtual Higgs fields. A simple ansatz for the field configuration describing the 
two bubble collision is 

{0 if v w t < z < —v w t, t < 

if -v w t < z < v w t, t > (13) 

v (T ) otherwise, 

where we have approximated the two colliding bubbles as plane-symmetric walls 

moving along the z direction with velocity v w . The Fourier transform of (13) 

yields 

«<*> = Pn^- (") 

kl - k 2 z vl 
From (12) we can estimate the number of particles x's produced per unit area 

/ d j^ lH{ko > h) ? Hkl - ® 9 $ ' kl - Am i) ■ (15) 



M r r dk z dkn 
~A~ J (2tt) 5 



The computation simplifies assuming the step-function configuration (13) and 
the scattering to be elastic [20]. These approximations give accurate results 
for k z < A^ 1 and k < (A^ 1 ^). Note that the slow power law-fall at large 
&Q — k^v^ in expression (14) is due to the fact that h(t, z) is discontinuous. For 
more realistic thick walls , Ti(ko, k z ) would cut off exponentially for k z > A" 1 and 
^^(A^^)" 1 [20]. 

If we consider a typical Yukawa coupling of the Higgs field to fermions x' s ) 
C = gh\X) from eq. (15) we obtain 

From the last expression we can argue that there is the possibility of producing 
particles with mass approximately up to 7A _1 . Hence, the value of 7 becomes 
crucial in our discussion. For our mechanism for AL production to work it 



is important to be able to generate particles with mass larger than rrih through 
bubble collisions and, thus, we need a large value of 7. The value of this parameter 
has been the focus of interest of a few recent analyses. In ref. [21] it is claimed 
that bubbles can expand with v w close to the speed of light. However Dine et al. 
[11] have identified relevant mechanisms wich may slow down v w . Their result 
depend on microsopic conditions for themalization processes in the vicinity of the 
wall and a more complete analysis of the Boltzmann equation to determine the 
departure from equilibrium in the particle phase-space density across the wall is 
needed [21]. Moreover, even a small gradient in temperature across the wall can 
be crucial for accelerating the wall itself, see Enqvist et al. [18] and ref. [22]. 

The estimate of fa comes from the direct comparison of M j A given in (16) 
and that calculated from (10) 

^~J^Af^V^fa~v{T ffa, (17) 

where the last equality is readily obtained using eqs. (4) and (9). We obtain 

We have now all the factors entering the expression (11) for the number density 
of particles x's. Then, AL produced in x decays is 

where 

r (x - k) - r (1 - u 



£ L = E B i 7u~~ 7n "» (2°) 

i r (x -»• all) 

Bi being the branching ratios of the decay mode x "~ *■ h an d 9*s (To) = 106.75 
counts for the relativistic degrees of freedom contributing to s, the entropy den- 
sity. In (19) we have taken into account that the "reheating" temperature is 
~ To. 



Using eq. (8) and taking a ~ 10 2 , u ~ 0.2, and At ~ 2 x 10 2 , we obtain 
from eq. (19) 

AL « 10- 2 £ L & 7" 1 . (21) 

As we previously mentioned, we consider i?-parity explicitly broken in the lep- 
tonic sector [10] as the concrete framework for the illustration of our ideas. The 
relevant L-violating terms in the superpotential are 

Wal^o = ~Kjk Li,Lj e c k + X ijk LiQjd c k , (22) 

where i,j, k are generation indices, L and Q are the left-handed lepton and 
quark doublet superfields, and e c and d c denote (the charge conjugate of) the 
right-handed lepton and charge —1/3 quark singlet superfields. 

There already exist several experimental constraints on the A and A' couplings 
in (22) [23] . The strongest bounds by far are those obtained from the cosmological 
considerations on the survival of the cosmic AB [24]. If we ask for the L-violating 
interactions due to the couplings (22) to be constantly out of equilibrium until 
the weak scale, so that they cannot wash out any (B — L) asymmetry previously 
generated, a severe upper limit A, A' ^ 10~ 7 , independently from the generation 
indices is obtained^. In our case these limits are no longer valid since we ask for 
the L-violating interactions in (22) to be still operative at the weak scale in 
order to give rise to the late AL. Other phenomenological and astophysical 
limits [26] exist, though they are substantially less tight, and we shall take them 
into account. 

The scenarios we have in mind for the generation of AL goes as follows. We 

assume the lightest supersymmetric particle (LSP) to be the lightest neutralino 

3 This is true only if one allows for the violation of all partial lepton number symmetries 
in (22). It was pointed out that one can evade these bounds if one of the lepton numbers 
L e , Lp, L T , say L e , is preserved in (22), so that iB — L e is not washed out [25]. 



X- Differently from what occurs in MSSM, the LSP can now decay through the 
^-violating couplings in (22). Since x possesses decay channels with different 
L- numbers, a net AL can develop in its decay if CP is violated and the out of 
equilibrium condition is implemented. CP violation can be easily provided by the 
complexity of the A and A' couplings in (22). This is a remarkable advantage with 
respect to late AB production in SM where CP violation provided by the CKM 
phase seems to be too small. Concerning the out of equilibrium requirement, this 
is automatically satisfied if the x's are produced in the bubble collision through 
their couplings to h. 

X is generally given by the superposition of the neutralinos of the model 
(notice, in particular, the presence of the singlet fermion N, in the case of one 
singlet superfield N) 

X = Z U W 3 + Z 12 B + Z^Hl + Z U H° 2 + Z l5 N. (23) 

As shown in ref. [27], when the LSP is x ~ \ \J^\ + -^2) — 775^' ^ can De ver y 
massive, m^ = 0(500) GeV. In this limit, the coupling of x to h reads 

C = Z 15 (Z 14 cos (3 + Z 13 sin (3) h^c—^X + h.c. (24) 

The coefficient g in eq. (18) is then 

g = V2Z 15 (Z u cos (3 + Z 13 sin (3) , (25) 

where y/2 takes into account that x is a Majorana fermion. 

The dominant decay channel of the LSP is expected to be that of fig. 1 0, 
whose thermally averaged partial width is 

<r (x - 1 k d%)) = r(x^tk d%) Kl { p^\ , (26) 

A 2 {m^/T ) 



4 If the i?-couplings are nonvanishing, LSP can decay also into two body final states via 
one-loop diagrams, for instance into WU, but these channels are suppressed by a factor O(f00) 
in comparison to that of fig. 1., see Campbell et al. [24]. 
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where K n denotes the modified Bessel functions and the prefactor in front of the 
usual (zero-temperature) decay rate can be interpreted as a time dilation factor. 
A straightforward calculation leads to 

where, in our case, Z u = |. 

The expansion rate of the Universe is given by H ~ 17T 2 /Mp, where Mp ~ 
1.2 x 10 19 GeV is the Planck mass. Neutralinos can decay after their production 
if, roughly speaking, T > H, which translates into the bound 

IV I > Q 2 x If)"* / 500GeV \ 5/2 / T X / m h \* ZlOOGeVx 



"x 



The interference of the diagrams of figs. 1 and 2 gives rise to a non vanishing El 
of order 

iv ( y* y y* y ) 

rvj Z-tikmst \/Htm/^ stm s3fc i3fc ) fOCA 

167r Eik\K 3 k\ 

Actually, this expression must be multiplied by a suppression factor which takes 
into account that part of the neutralinos produced out of equilibrium by our 
mechanism may thermalize before decaying, e.g. through the xt — > xt scattering, 
and, consequently, they cannot give contribution to El- A rough estimate for such 
a suppression factor is given by (r (x — »■ t U d c k )) / (r (t X ~ * t x)) ~ 10 3 | A' 13fc | . 

In order to prevent a washing out of the lepton asymmetry, we impose that 
AT = 1 scatterings induced by A and/or A' couplings are out of equilibrium at 
To. The predominant process is shown in fig. 3. Since To < m^, this process has 
a thermally averaged rate 
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r AL=1 = n EQ (T)(av) ~ 32^3)1^3,1 -^-<xw x 

x (^l) 2 - * e -«*/*b, (30) 

where £(3) ~ 1.202 is the Riemann zeta function and n^g denotes the equilibrium 

number density. Requiring Tal=i ^ # imposes 

iw i^xo , ft -4/ m <L \ 2 /500GeV\ /150GeV\ 1/2 /100GeV\ 

I A^l £ 4.2 x 10 (^) ^__ j (^_) (__) . (31) 

A similar upper bound is found requiring that inverse decays are out of equilib- 
rium. The products of neutralino decays will be brought to thermal equilibrium 
at rates of order To. At this temperature, anomalous (B + T)-violating processes 

4n v(T ) 

are in equilibrium, since T sp h a ieron — ®wTe B w T o > H for At ~ 2 x 10 and 
the lepton asymmetry is converted into a baryon asymmetry AB = Ub/s, whose 
present value is {T tod ~ 2.7°K) [28] 

AB ^ _|1^1M AT, (32) 

61 g*s{To) 

where g* s (T to d) / g* s (T ) takes into account the increase of entropy since the 
EWPT up to now. 

From eqs. (21), (29) and (32), we argue that a baryon asymmetry 4 x 10 -11 ^ 
AB < 5.7 x 10~ n [29] can be produced for resonable values of the Tt-violating 
couplings. For instance, taking mf L ~ 5 TeV, 7 ~ 10 2 , At ~ 2 x 10~ 2 and all the 
A' couplings of the same order, we obtain A'~8x 10~ 3 . Note that this value is 
in agreement with the upper bounds given in [26]. 

Notice that our scheme makes use of the Tt-parity violating couplings in the 
leptonic sector. Previous proposals to use .R-parity violation for late AB pro- 
duction invoked Tt-breaking in the baryonic sector within contexts quite different 
to implement the out of equilibrium condition [30] . 

12 



In conclusion, we have proposed a scheme for late baryon asymmetry in a first 
order EWPT through the interplay of lepton violating interactions and anomalous 
B- and L-violating quantum effects: i) AL is generated in the decay of the 
LSP produced out of equilibrium in bubble collision and ii) AL is converted 
into AB by sphalerons still operative after the EWPT. Obviously, a thorough 
discussion of several open issues of EWPT is eagerly awaited for to provide a 
more comprehensive and detailed scheme for late AB production along the lines 
suggested in this paper. 
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Figure Captions 



Figure 1 . The tree-level Feynman diagram for the decay x ~~ *■ thd% induced by 
the \' i3k couplings. 

Figure 2. The one-loop Feynman diagram for the decay x ~~ * t k d c k , whose inter- 
ference with that of Figure 1 gives rise to a non vanishing lepton asymmetry. 
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Figure 3. The Feynman diagram for the AL = 1 scattering li d c k — > tx, which 

can wash out any lepton asymmetry. 



